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I .     INTRODUCTION 

Free  Electron  Lasers  (FEL)  may  have  an  impact  on  the 
United  States  Navy,  not  just  as  a  weapon,  but  in  many  fields 
of  science  such  as  Medicine  and  Physics.  The  tunability  and 
high  power  output  of  a  FEL  makes  it  useful  for  a  broad  range 
of  applications. 

As  a  weapon,  the  tunability  allows  the  user  to  select 
the  frequency  that  will  propagate  most  efficiently  through 
the  atmosphere.  Also,  the  FEL  has  the  potential  to  produce 
the  power  level  needed  for  use  as  a  weapon. 

The  tunability  allows  a  doctor  to  select  not  only  the 
wavelength,  but  also  the  pulse-repetition  rate  of  the  laser 
during  surgery.  The  FEL  is  able  to  take  advantage  of 
resonant  modes  in  a  molecule  to  allow  the  doctor  to 
precisely  remove  millimeter  layers  of  skin,  fat,  muscle, 
cornea,  or  tumor  with  no  collateral  damage  to  surrounding 
tissue.  The  FEL  "scalpel"  may  revolutionize  the  operating 
room. 

The  FEL  is  also  the  most  viable  path  for  achieving  a 
coherent  x-ray  source.  There  are  several  projects  working 
towards  a  high  gain  FEL  that  can  achieve  1  angstrom  x-rays. 
The  uses  for  coherent  x-ray  sources  are  numerous.  One  use, 
which  has  a  direct  effect  on  the  Navy,  is  to  probe  the 
reactions  that  occur  in  the  initial  moment  of  a  nuclear 
blast . 


There  are  numerous  other  uses  for  the  FEL.  Almost  every 
field  of  science  is  finding  new  applications  for  the  FEL, 
from  probing  light  interactions  with  biological  functions  to 
studying  Quantum  Electro-dynamics. 

This  thesis  covers  several  different  FELs  which  have  a 
wide  range  of  applications.  The  second  chapter  of  this 
thesis  covers  the  background  theory  that  will  be  used  in  the 
rest  of  the  thesis [reference  (1)]  .  In  Chapter  III  I 
calculate  analytically  the  height  of  a  tapered  separatrix. 
In  Chapter  IV  I  analyze  various  high  gain  FEL ' s  and  use  some 
dimensionless  parameters  to  develop  a  simple  model  for  each. 
The  final  chapter  I  simulate  the  Darmstadt  FEL  using  a 
computer  model  and  analyse  the  results .  This  analysis 
consists  of  simulations  at  various  desynchronisms  that 
determine  power  and  gain  expected  in  the  FEL.  There  are 
instabilities  that  may  be  experimentally  observed  as  well. 


II.   BACKGROUND  THEORY 

The  physics  behind  the  Free  Electron  Laser  (FEL)  is 
based  upon  both  spontaneous  and  stimulated  emission  from  an 
electron.  The  two  major  components  of  the  FEL  are  the 
relativistic  electron  beam  and  a  static  alternating  magnetic 
field.  The  relativistic  beam  of  electrons  is  accelerated  in 
the  transverse  direction  by  the  magnetic  field.  This 
acceleration  initially  causes  spontaneous  radiation  from  the 
electrons.  Due  to  relativistic  effects,  the  emitted  light  is 
radiated  into  a  cone  along  the  axis  of  the  laser.  The  key  to 
an  operational  FEL  is  the  interaction  between  the  electrons 
and  the  light  field.  The  kinetic  energy  of  the  electron  beam 
is  converted  to  the  radiation  energy  of  the  optical  field. 
The  electron  pendulum  equation  and  optical  field  equation 
describes  the  coupling. 

A.     INITIAL  CONDITIONS 

A  classical  approach  has  proven  to  be  the  most  fruitful 
and  experimental  results  agree  with  the  classical 
results,  [reference  1] 

1.    Undulator  Field 

The  magnetic  field  of  the  FEL  is  created  by  a  chain  of 
fixed  magnets,  called  an  undulator.  There  are  two  basic 
types  of  undulators :  helical  and  linear.  The  linear 
undulator  is  a  series  of  bar  magnets  with  the  polarization 


in  alternating  directions,  so  that  an  alternating  magnetic 
field  is  set  up.  Figure  1  depicts  a  linear  undulator  with 
the  z  direction  is  along  the  axis  of  the  laser.  The  magnetic 
field  for  the  linear  undulator  is  given  by 

B  =  [bx,  By,  B.J  =  B_[o,  sin(kcz)cosh(k0y),  cos(k0z)  sinh(k0y)] ,     (II-l) 


where  k0  =  2tt  /  X0    is  the  undulator  wavenumber 
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Figure  II-l  Linear  Undulator  (after  reference  7) 

The   helical   undulator   consists   of   two   concentric 

currents  that  create  a  magnetic  field  that  rotates  spatially 
along  the  length  of  the  undulator.  For  the  helical  undulator 
the  magnetic  field  is  given  by 

B  =  Bo[cos(k0z),sin(k0z),0]  (II -2) 


2.    Initial  Optical  Field 

To  begin  analyzing  FEL  operation  we  need  to  assume  that 
the  FEL  is  already  operating  and  there  is  already  a 
radiation  field  present.  The  optical  field  for  the  linear 
undulator  is  assumed  to  be  a  plainer,  plane  parabolic  wave 
with    the    electric    and   magnetic    fields    given   by 

Es  =  E[cos  i|/,0,0]   and   Bs  =  e[o,  cos  i|/,0]   respectively,   where 

i|/  =  kz  —  cot  +  <j>  .  The  electric  field  amplitude  E  is  given  in 

cgs  units  with  (j)  being  the  optical  phase  and  k  =  2%  /  k     is 

the  optical  wave  number,  co  =  2nc  /  A,  is  the  optical  angular 
frequency. 

The  initial  optical  field  for  the  helical  undulator  is 
assumed  to  be  a  circularly  polarized  plane-wave.  Hence  the 
electric     and     magnetic     fields     are     given     by 

Es  =  E[cos  i|/,- sin  i|/,0]  and  Bs  =  E[sin  vy,  cos  i|/,0] ,  respectively. 

B.  ELECTRON  DYNAMICS 

1.    The  Electron  Phase 

Before  moving  on,  we  need  to  define  the  electron  phase. 
The  electron  phase  is  the  key  to  coupling  the  optical  field 
with  the  bunching  of  the  electron  beam.  In  the  combined 
undulator  and  light  field  the  electron  phase  is  defined  as 

C,  (t)  =  (k+k0)  z  (t)  +G)t .  The  electron  phase  definition  involves 
two  fixed  wave  numbers  k  =  2n  I  X  ,  describing  the  light  wave 


and  k0  =  2tt  /  Xc  describing  the  undulator,  with  the  only 
variable  being  z(t) .  The  electron  phase  follows  the  electron 
position.  The  time  spent  by  the  electron  in  the  undulator  is 

At=L/P0c^L/c~10"8s,  where  L  is  the  length  of  the  undulator 

and  the  speed  of  the  electron  (v2=cPz«c)  is  approximately  the 

speed  of  light.  This  time  is  extremely  small,   so  it  is 

easier   to   introduce   a   dimensionless  time  x       and   a 

dimensionless  time  derivative 

cB_t    ct    ,  /  °  \    d(.  .) 
t  =  -^-  «  —  and  (.  .)  =  -^.  (II-3) 

L      L       v  J  dx 

The  dimensionless  time  x  ranges  from  0  at  the  beginning  of 

the  undulator  to  1  at  the  end  of  the  undulator  (0<t<1) .  The 

electron  phase  velocity  is  then  defined  as  the  derivative  of 
the  electron  phase 

v(x)  =  ^Q  =  i  =   L[(k  +  k0)32  -  k]  (II-4) 

dx 

To  study  the  operation  of  the  FEL,  it  is  most 
informative  to  plot  the  phase  velocity  versus  the  electron 
phase  in  what  is  called  phase-space.  This  will  be  discussed 
further  in  subsequent  parts  of  this  chapter. 

2 .    Wavelength  of  Light  in  FEL 

The  traditional  laser  requires  a  bound  electron  to 
transition  to  a  lower  energy  level  and  radiate.  In  the  FEL, 
a  relativistic  beam  of  free  electrons  amplifies  the  light 


wave.  The  undulator  wavelength  ^0  is  Lorentz  contracted  to 

X0     =   X0  /  y  in  the  electron  frame.  The  emitted  optical  field 
is  Doppler  shifted  is  also 

X   =  (l  -  P2)yA.  *  2yX  .  (II-5) 

At   resonance,   there   is  optimal   energy   exchange.   When 

X0     =   X'     in  the  beam  frame.  When  converted  to.,  the  lab,  frame 

the  resonance  condition  becomes 

»      1  +  K2 

A,  =  XQ  —  (II-6) 

2y 

where  the  dimensionless  undulator  parameter  K  is  defined  as 

eBX0      el 

K  =  °—   =  -,  (II-7) 

27i  mc     k0mc" 

B  is  the  RMS  magnetic  field,  e  is  the  magnitude  of  the 
electron  charge,  m  is  the  mass  of  the  electron  and  c  is  the 
speed  of  light. 

3 .    Electron  Paths 

To  analyze  the  electron  trajectories  in  the  FEL,  we 
will  look  at  the  linear  undulator.  The  magnets  are  aligned 
along  the  y  direction.  An  electron  that  moves  off  axis  in 
the  y-direction  feels  a  stronger  magnetic  field  because  it 
moves  closer  to  a  magnet.  The  electron  is  pushed  back  toward 
the  axis  creating  a  natural  focusing.  The  focusing  of  the 
undulator  causes  the  electrons  to  wiggle  over  many  undulator 
periods  in  the  focused  direction.  The  wiggle  is  called  a 


betatron  oscillation.  To  analyze  betatron  effects,  we  will 
consider  a  perfectly  injected  electron  in  a  linear  undulator 
that  focuses  in  the  y  direction.  Perfect  injection  means 
that  the  electron  is  injected  into  a  perfect  sinusoidal 
trajectory  with  no  angular  spread. 

The  three  spatial  components  are  of  the  Lorentz  force 
equation  with  the  magnetic  field  for  linear  undulator  from 
equation  (II-l)  as  follows 

p.  = [p.  sin(k0z)cosh(k0y)  -  B  cos(k0z)  sinh(k0y)J ,        (II-8) 

ymc 


P„  =  [p..  cos(krz)sinh(kcy)],  (II-9 

ymc 


p,  =  - -^— [P.,  sin(k0z)cosh(k0y)]  .  (11-10) 

ymc 

In  the  next  section,  we  will  derive  these  equations  in 

greater  detail  for  the  helical  undulator.  The  Px  term  has  an 

exact  solution,  and  perfect  injection  implies  the  constants 
of  integration  are  set  equal  to  zero,  so 

pv  =  cos(k0z)  cosh(k0y) .  (11-11) 

Y 

Substituting  Px  into  the  other  two  equations  leads  to 

cK2k 
(3   = cos2(k0z)sinh(2k.y),  (11-12) 

Y~ 
and 

cK2k 
pz  =  —  sin(2k0z)cosh2(k0y).  (11-13) 


Since  the  y  and  z  components  are  coupled,  we  add  the 
restriction  of  perfect  injection  in  the  y  direction  and 
analyze  the  motion  in  the  xz-plan.  Hence, 

a  V2K     /n   >. 

P,  = cos(k0z)  (11-14) 

Y 

Py(0)  =  y(0)  =  0  (11-15) 

P,  =  r1  sin(2k0z)  (11-16) 

Y~ 

In  the  typical  operating  FEL  K«l  and  y>>l,  so  that 
K/y<<l.  Expand  in  powers  of  (K/y)  ,  so  the  x  and  z  components 
become 

z(t)  =  cp_t sin(2k0cp,t)  +  ...,  (11-17) 

4k0y2    v     ' 

x(t)  = sin(k0cp,t)  +  .  .  .  ,  (11-18) 

where  Pz  is  the  average  velocity  of  the  electron  beam  in  the 
undulator.  Then  converting  to  dimensionless  time  x  gives 

—  K2X 

z(x)  =    p,Li 2.  sin(47iNx)  +  .  .  .  ,  (11-19) 

87ry" 

KX 

x(t)   =    -  -=$-  sin(2;rNT)  +  .  .  .      .  (11-20) 

V27iy 

From  these  two  equations,  we  can  see  that  the  amplitude  of 
the  transverse  oscillation  is  KA,0  /  V27ry  where  the  electron 


oscillates  once  every  undulator  period.  In  the  z  direction, 
the  electron  oscillates  twice  every  undulator  period. 
Having  looked  at  the  oscillations  in  the  x  and  z  directions, 
we  now  go  back  and  assume  non-perfect  injection  to  see  what 
happens  in  the  y  direction.  The  oscillations  in  the  y 
direction  are  called  betatron  oscillations.  Assuming  small 
oscillations   so   that   k0y<<l   and   sinh(2kcy)  «  2k0y   we.  can 

rewrite   equation   (11-12)   and   convert   to   dimensionless 
parameters  to  get 

2K2L2k  2 
y(x)  = r-5-  cos2(27iNT)y(x)  (11-21) 

Y" 
Averaging  over  several  undulator  periods,  since  one  betatron 
oscillation  is  typically  over  many  undulator  periods,  we 
have 

K2L2k  2 

y(x)  = r^-yW  =  -oVyd:).  (11-22) 

Y 

Equation   (11-22)   is  the  equation  of  motion  for  a  simple 

harmonic   oscillator.   The   solution   is   a   well -studied 

sinusoidal  wave,  with  the  dimensionless  betatron  frequency 

cop=KLk0/y=2TcNK/y.  The  number  of  betatron  oscillations  along 

the  undulator  is 

NK 
NB  =  .  (11-23) 

Y 


10 


A  typical  value  is  Np«l ,  since  N-100,  K=l,  and  y=100  are 
typical  values.  FELs  with  larger  values  will  be  discussed 
later . 

C.     ELECTRON  PENDULUM  EQUATION 

The  dynamics  of  the  FEL  may  be  described  very  simply  by 
the  electron  pendulum  equation.  The  first  step  in  deriving 
this  equation  is  to  consider  an  operating  FEL  where  an 
optical  field  is  already  present.  This  means  the  electron 
experience  forces  from  both  the  optical  field  and  the 
undulator  field. 

The  helical  field  is  much  simpler  to  work  with 
mathematically  and  the  final  formulas  can  be  converted 
easily  to  the  linear  case.  Hence,  we  will  just  use  the 
helical  magnetic  field. 

The  forces  from  the  undulator  field  are  governed  by  the 
Lorentz  force  equations  (cgs  units) , 

di  =  d(ymcp]  =  _e(g  +  pxi)^  (I 

dt     dt       v       ; 

and  the  energy  equation 

dlymc2 )  -  , 

-^ '-  =  -e   cB-E,  (11-25) 

dt 

where 

^r    =   1  -  B2.  (11-26) 
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The  momentum  is  p=ymv=ym|3c,  its  energy  is  ymc2,  and  the 
charge  magnitude  e  =  |e|.  We  have  five  equations  and  four 
unknowns  x,  y,  z  and  y. 

Substituting  the  magnetic  field  from  equation  (II-2) 
into  equations  (11-24)  ,  (11-25)  and  (11-26)  yields  the 
following  equations 

Md  =  _  _•.  [E(1  _  P:Xcos  Vi_  sin  v,0)  + 

dt  mc  /  [ll-z  / ) 


p2B(-  sin(k0z),  cos(k0z),0)] 
= [e((3.,  cos  \\j  -  (3y  sin  \\i)  + 


d(y(U 

dt  mc  L  v  '" 

b((3:.;  sin(k0z)  -  py  cos(k0z))] 


— —   = e[Px  cos  \\j  -  py  sin  \\f] , 


11-29 


dt      mc 
where   (3X  =  (3,.  +  |3y   and   vj;  =  kz  -  cot  +  <|)  .   There   are   more 

equations  than  unknowns,  so  we  can  ignore  the  P2  component. 

The  electrons  are  traveling  near  the  speed  of  light 
along  the  z  axis,  so  we  can  estimate  that  Pz«l .  This  means 

that  E(l  -  P7)  «  BP_  ,  and  the  first  term  of  the  transverse 
components  of   (11-22)  can  be  ignored.  Thus, 


[PZB(-  sin(k0z),  cos(k0z),0)],  (11-30) 


dt       mc 

Which  implies  that  the  static  undulator  field  controls  the 
transverse  motion  of  the  electron  almost  entirely.      We 
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can  integrate  (11-30)  by  inspection,  yielding  the  following 
results 

Pi  = [cos(k0z),  sin(k0z),o]  .  (11-31) 

y 

Substituting  (11-31)  into  (11-29)  yields 

eKE    f„        ,x 

y  =  cos(<;  +  ({))  (11-32) 

ymc 

where    C,  +  <j>  =  \|/  +  k0z  =   (k  +  k0)z  -  cot  +  <j>  .    From    (11-31) 

K2 
Px2    =   —       ,  (11-33) 

y 

and   substituting    into    to    (11-26) 

\    =    1   -   (p)2    =    1   "   Pz2    "  P/    =    1   "   P,2    "  \  (11-34) 

Y  Y 

so    that 

(l  +  K2)^r    =    1  -  P,2  .  (11-35) 

y^ 

Taking    the    time    derivative    of     (11-35), 

Y  Y2P  3 

I   =     r  PzPz?  (11-36) 

y         1  +  K" 

From  the  derivative  of  the  electron  phase,  we  have 

4  =  (k  +  k0)z  -  co  =  (k  +  k0)cP  -  co  ,  (11-37) 

and 

I   =   (k  +  k0)cpz,  (11-38) 

so  that 
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(k  +  k0)c  ' 
Substituting    into     (11-36),    we    get 

y   y2PA      y2P^ 

-  = -  i n? \~  •  H-40 

y         1  +  K"         (l  +  K~Xk  +  kc)c 
Now   equating    (11-32)    and    (11-40),    we   get 

7 j£r* v-   =   — —  cos  (£  +  <j>)  ,  (11-41) 

(1  +  K~  )(k  +  k0jc  mc 

and  solving  for  C,  ,  we  have 

■■         (l  +  K2)  (k  +  k0)c  eE0K        ^ 

C,    =  A : r— —  COS(C  +  (j)).  (11-42) 

y  p:  y"mc 

Since  k>>k0i  we  can  use  the  approximation  (k+k0)~k.  From 

equation  (II -6)  co0  =  ck0  =  ck(l  +  K2)/2y2,  and  using  (3Z«1, 

equation  (11-42)  reduces  to 

<;  =  \ cos(C  +  <|»  (11-43) 

y^mc 

Converting  to  dimensionless  time  with  dt=Ldx/c,  and  defining 

the  dimensionless  optical  field  magnitude  as 

■  ,    47iNeKLE 

a  =  — - — — ,  (11-44) 

y^mc" 

The  pendulum  equation  becomes 

o  o 

<;  =  |a|  cos  (<;  +  (j))  .  (11-45) 

The  pendulum  equation  (11-45)  is  the  same  form  for  a  linear 
undulator.  However,  for  the  linear  undulator,  the  optical 
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field  |a|  is  modified  by  substituting  K [J0  (^)  - Jx  (<*;)  ]  for  K, 
where  ^=K2/2  (1+K2)  and  J0  (£,)  and  J1  (£,)  are  Bessel  functions. 

D.     OPTICAL  WAVE  EQUATION 

The  pendulum  equation  describes  how  individual 
electrons  are  affected  by  the  laser  field.  To  describe  how 
the  laser  field  is  affected  by  the  beam  of  electrons,  we 
next  derive  the  optical  wave  equation.  Begin  with  the 
inhomogeneous  wave  equation  for  the  optical  vector  potential 

A, 


V2-^ 


c~  at' 

where 


471  - 

A  = J±l  (11-46 

c 


Ee    = =    E[cos  \\if-  sin  \\/,0] ,  (11-47) 

c  dt 

Bs  =  V  x  A  =  [sin  i|/,  cos  Vj/,0] ,  (11-48) 

\\j   =   kz  -  cot  +  (j)  .  When  the  amplitude  and  phase  of  the  laser 

beam  change  slowly  with  time,   we  can  write  the  vector 
potential  in  the  form 

A(z,  t)  =  —  E(z,  t)  [cos  (V|/)  ,-  sin  (\j/)  ,0]  .  (11-49) 

CO 

The  form  in  equation  (11-49)  assumes  that  the  electron  beam 
is  small  compared  to  the  waist  of  the  optical  mode,  so  that 

the  transverse  direction  is  not  important.  A  is  dependent 
only  on  z  and  t.  The  spatial  derivative  is 


15 


d2A 

dz2 


1  dE 
k  dz 
1   d2E 


dz 


[cos  (v|/)  ,-  sin  (\|/)  ,0]  + 


k 

dz2 

1 

dE 

k 

dz 

E 

d2§ 

[sin  (i|/)  ,  cos  (vj/)  ,0] 


+ 


k  + 


5z 


[cos  (\j/)  ,-  sin  (v|/)  ,0]  + 


(II-50 


k  dz' 

E 


[cos  (v|/)  ,-  sin  (if/)  ,0]  + 


a^ 

dz 


[-  sin  (\|/)  ,-  cos  (v|/)  ,0] 


and  the  time  derivative  is 


co  5ZA         dE 


c  dt~         dt 


~dt 


co 


a2E 


+ 


5E 

at 


at 


-    CO 


[cos  (v|/)  ,-  sin  (vj/)  ,0]  +  — -  [sin  (\|/)  ,  cos  (\j/)  ,0] 

at" 

a<|>2 


[cos  (\j/)  ,-  sin  (v|/)  ,0]  +  E  — *-r  [cos  (\|/)  ,-  sin  (i|/)  ,0]     (11-51' 

at1 


+  E 


a^_ 

at 


CO 


[-  sin  (v|/)  f-  cos  (\|/)  ,0] 


From     the     assumption     that     phase     and     amplitude     varies 
slowly,      dEjdz     <<     kE  ,      d§/dz     <<     k§  ,      dE/dt     <<     coE  ,     and     dfy/dt 

<<co(j)  .  From  these  approximations,  we  find  that 


v*  -  i  5 


c2  at: 


A*2 


dE     J_aE 

az    c  at 


[cos(v|;),sin(v|/),0] 


11-52 


+  2E 


a<j)   j_a^ 
az    c  at 


[-sin(v|/),-cos(i}/),0]  = Jx 

c 


To  further  simplify,  we  use  a  coordinate  transformation, 
where  the  position  z  moves  along  with  the  optical  field  at 
the  speed  of  light  c,  so  that  z  =  z  +  ct  .  Also,  we  convert 
to  dimensionless  time,  x=ct/L  to  get 
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d_ 
dz 


d 


and 


]_dE 
L  dx 


cdt)        cdt    Ldx 


[cos(v|;),sin(v|/),o]+  2E 


(11-53) 


L  dx 


[-  sin(i|/),-  cos(i|/),Oj  t 


471- 

J 


11-54 


Each  electron  adds  to  the  total  current,   so  the  current 
density  contribution  from  the  ich  electron  is  summed  to  get 


x  -r, 


11-55) 


:(3) 


where  the  r,  is  the  trajectory  of  the  i  electron  and  8 
is  a  three  dimensional  Dirac  Delta  Function.  Substituting, 
the  transverse  velocity  p±  =  -K  /  y(cos(k0z),  sin(k0z),0)  into  the 
current  gives 


-  —  Jj.~y! [cos  (koz)  i  s^n  (k0z)J  5(3)(x  -  ri ) 

c     i     Y 

Therefore,    the   wave   equation   is 


11-56 


dE 
dx 


[cos  (i|/)  ,  sin  (\\f)  ,0]  +  E 


dx 


[-  sin  (\\i)   -  cos  (v|/)  ,0]  = 


X"  27teKL[cos  (k0z)  ,  sin  (k0z)]  8(3)(x  -  rj 

i    y 


(11-57) 


Both  sides  of  (11-57)  averaged  over  a  volume  element  dVc, 
which  is  smaller  than  the  coherence  volume,  yet  larger  than 
a  wavelength  of  the  light.  The  current  density  over  this 
small  volume  element  is  constant,  so  the  summation  over  i 

electrons   V  (-  27teKL  /  y)[cos  (k0z)  ,  sin  (k0z)J  8(3)  can  be  replaced 
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with   a   charge   density  p   times   the   average   of   sample 

electrons  in  the  beam,  the  average  is  symbolized  by  (  )  .  So 
now  we  have  two  uncoupled  equations 

dE  27iepLK 

dx  y 

d§         27repLK 


(cos  (<;  +  (j)))  (11-58) 

(sin(<;  +  (j)))  (11-59) 


dx  y 

where    ij/+k0z=(^+<j).    Defining   the    dimensionless    current 

j  =  87re"pNK2l7/yJm  c"  ,  and  recalling  the  dimensionless  optical 

field   |a|  =  47iNeKLE/Y2mc  ,   the   change   in   electric   field 
amplitude  becomes 

-^  =  a  =  -j(cos((^  +  <}>))  .  (11-60) 

Since  the  electric  field  is  complex,  then  a  =  |a|e1<J>  ,  and  we 
arrive  at  the  optical  wave  equation 

a  =  -j(e"ic)  .  (11-61) 

Equation  (11-61)  is  derived  for  a  helical  undulator. 
However  the  same  form  of  the  equation  applies  to  a  linear 

undulator  provided  K  is  replaced  by  K  [  J0  (£)  -  J1  (£,)  ]  . 

E.     GAIN  MECHANISMS 

For  the  FEL  to  be  operational,  the  there  must  be  an 
increase  in  power  of  the  optical  field  each  pass  through  the 
undulator.  There  are  two  regions  of  optical  gain:  the  low 
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gain  region   (j<7r)   and  the   high  gain  region   (j>>7t).   To 
understand  the  mechanisms  behind  gain  we  first  develop  the 
concept  of  an  electron  phase  space  plot. 
1 .    Phase  Space 

The  electron  phase  space  is  a  plot  of  an  electron's 
phase  velocity,   v(i)  =  t,  =  h[(k   +  k0)3z  -  k]  and  the  electron 

phase  C,  (t)  =  (k+k0)  z  (t)  +0)t .  In  the  low  gain  case  a  *  0,  which 
means  that  equation  (11-61)  evolves  slowly  so  that  |a|  and  <j) 
are  approximately  constant.  Multiplying  both  sides  of  the 

o 

electron  pendulum  equation  (11-45)  by  C,    yields 

£  £  =  |a|  cos  (<;  +  (j>)  t,   =   v  v  (11-62) 

or 

v  v  -  |a|  cos  (£  +  <j>)  4  =  0     ,  (11-63) 

so  that 


d 

dx 


—  ~  |a|  sin(^  +  <j>) 


=  0 .  (11-64) 


This  means  that  v2  -  2|a|  sin((^  +  ()))  =  cons  tan  t  ,  and 

v2   =  v02  -  2|a|[sin(^  +  <\>)  -  sinfe)].  (11-65) 

Where  v0=v(0)  and  ^0=^(0)  are  the  initial  phase  velocity  and 
phases.  Equation  (11-65)  is  a  statement  of  energy 
conservation  in  the  dimensionless  coordinates  C,    and  v.  The 
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equation  also  identifies  paths  in  the  phase  space,  v(^), 
that  describes  the  electron  motion.  In  Figure  II -2,  several 
phase  space  paths  are  plotted.  8  refers  to  taper  and  will  be 
discussed  in  the  next  section. 


FEL  PHASE  SPACE  PATHS 
a=l         (j)=0         8=0 


-7C/2 


c 


37U/2 


Figure  II -2  Phase  Space  Diagram 
2 .    Separatrix 

Electrons  stay  on  v  {C,)     paths  in  the  C,-v    plane.  These 

paths  are  based  upon  the  initial  phase  velocity  v0  and 

initial  phase  <^0  of  the  electron  and  the  pendulum  equation. 
This  same  equation  describes  an  undamped,  undriven  pendulum 
where  C,    is  the  tilt  angle  and  v  is  the  angular  velocity.  The 
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critical  points  for  these  paths  occur  at  v=0  and  ^=7t/2±7in 
where  n=l,2,3.  The  stable  fixed  points  occur  when  Cs=n/2±2nn. 
In  Figure  II-2  the  stable  fixed  point  is  labeled  with  an 
"S"  .  The  unstable  fixed  points  are  at  v=0  and  C,=3iz/2±2nn .  In 
Figure  II-2,  the  unstable  fixed  points  are  labeled  with  a 
"U" .  At  these  fixed  points,  the  electron  remains  stationary 

in  phase  space.  At  (v,^)  =  (n/2, 0)  ,  the  pendulum  is  hanging  at 

the  bottom  of  the  arc.  When  |v|>2|a|1/2,  the  pendulum  has 
enough  angular  momentum  to  swing  continually  over  the  top 
around  in  a  circle.  In  phase  space,  we  can  draw  a  separatrix 
which  separates  the  lower  closed  orbits  from  the  over-the- 
top  open  orbits.  The  equation  for  the  separatrix  is 

vs2  =  2|a|[l  -  sin(<;s  +<())].  (11-66) 

In  Figure  II-2,  the  separatrix  separates  the  open  orbits, 
shown  as  a  gray  area,  from  the  closed  orbits,  shown  as  a 
white  area. 

3 .    Bunching 

Consider  a  collection  of  electrons  that  are  initially 
uniformly  distributed  in  phase  evolving  along  the  undulator. 

The  pendulum  equation  shows  that  for  -n/2  <  C,  <  n/2,  C,>0  and 
the  electron  phase  accelerates,  while  for  n/2     <    £,    <    37r/2, 

<^<0   and  the  phase  decelerates.   The  result   is  that  the 
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electrons  tend  to  bunch  about  ti/2  as  they  evolve  along  the 
undulator . 

4.  Changes  in  Phase  Velocity 

From  the  definition  of  phase  velocity 
v(i)  =  L[(k  +  k0pz  -  k],  we  find  the  change  in  beam  energy  due 
to  a  change  in  phase  velocity.  Taking  the  derivative  of  both 
sides  near  resonance  where  v=0,  k>>k0  and  N=h/X0,  we  get 
Sv  w  Lk6(3z  «  N?t0k8|32  .  From  equation  (II-6)  and  the  resonance 
condition,  X  =   A0(l  +  K2)/2y2    ,  we  find  that 

8v  «  47iN  —  «  2ttN  —  .  (11-67) 

Y        A. 

Thus,   an  increase  in  v  results   in  an  increase  of  the 

electron  energy,  ymc2,  and  the  optical  wavelength  X. 

5 .  Gain  Formula 

Gain  is  defined  as  AP/P,  where  AP  is  the  change  in  optical 

power  in  one  pass  through  the  undulator,  and  P  is  the  power 
at  the  beginning  of  the  undulator.  From  (11-67)  the  average 
change  in  energy  for  an  electron  in  the  undulator  is 

Aymc:  w  ymc2((v)  -  v0)/4ttN,  (11-68) 

where  (v)  is  the  average  final  energy.  The  number  of 
electrons  in  a  volume  dV  of  the  optical  wave  is  given  by 
pFdV  where  p  is  the  electron  density  and  F  is  the  filling 
factor.  The  filling  factor   F  is  the  cross-sectional  area  of 

22 


the  electron  beam  divided  by  the  cross-sectional  area  of  the 
light  beam.  The  radiation  energy  in  a  volume  dV  is  given  by 

2E2dV  /  8tc  in  cgs  units.  Invoking  the  definition  of  the 
dimensionless  field  a  and  the  dimensionless  current  j  ,  we 
arrive  at 

[PFdv][ymc2((v)  -  v0)/4ttn]        2Fj  ,  ,  (n    „. 

b    =    1 — ; ■ 1 -   — -  (v0    -  V).  (11-69) 

[2E^dV  /  8ttJ  a02   v  ' 

6 .         Low  Gain 

For  low  gain,  the  optical  field  changes  vary  slowly  so 

we  can  estimate  that  a  «  0  with  the  initial  conditions 
v(0)=v0/  a(0)=a0,  and  <^(0)=^0.  From  perturbation  theory,  we 
find  that 


a 
v  =  v 

v      v0 

V 


+  ^\  [sinfe0   +  v0t)  -  sinfc)]  + 


2 

v„3 


—  cos(2^0  +  2v0x)  -  cos(2^0)  +  cos(2v0t)  - 

4 


(11-70) 


1  -  v0x  sinfcjcosfc  +  v0x) 
Taking   take    the    time   average   of    the   phase   velocity,    we    find 

(v)   =   v0  +  ^2_  [2  cos(v0x)  -  2  +  v0t  sin(v0x)]  +  .  .  .       .  (11-71) 

vo 

The   gain   is    therefore 

G   =   -^f  (v0   -  v)   =   ^j  [2  -  2  cos(v0t)  -  v0t  sin(v0x)] .  (11-72) 


a0  vo 
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7 .    High  Gain 

In  the  high  gain  regime  the  current  is  very  large, 
j>>7i.  Unlike  the  low  gain  regime,  we  cannot  ignore  the  rapid 
evolution  of  the  optical  field  over  a  single  pass.  Combining 
the  expanded  pendulum  and  optical  equation  and  integrating 

over  all  the  initial  electron  phases,  jd^0  ,    we  get 

a(x)  =  ^  [dxVF(T')e-lu»Ta(i  -  x')  (11-73) 

where  F(t')=  jdqf(q)e~'qT  is  the  characteristic  function  of  the 
distribution  f  (q)  .  f  (q)  is  the  distribution  of  the  initial 
electron  phases  vi=v0+q  about  v0  and  Jdqf (q)  =  1  . 

The  low  current  gain  formula  can  be  obtained  from  this 

equation  by  considering  a  perfect  beam  where  f(q)=8(q)  and 

a(x  -  x')  ~  a0  ,   in  which  case  a0  can  be  moved  outside  the 

integrand.  Further  integration  results  in  the  low  gain 
formula  shown  earlier. 

To  find  the  solution  to  equation  (11-73),  we  need  to 

consider  a  perfect  beam  with  F(x')  =  l  and  a  resonant  optical 
field  v0=0.  Then,  taking  two  derivatives,  we  get 

a(x)  =  —  a(x).  (11-74) 

It  can  be  shown  that  the  solution  to  (11-74)  is 
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a(T)  =   £o.[e(3./2)-3(i+V3)t/2  +  eCj/2j">a-V3>/2  +  elJ/2j"',]  (II-75) 

Since  for  high  gain  j>>7i  the  dominate  exponential  terms  are 

the  fastest  growing  terms.  The  following  approximations  can 
be  made  for  the  optical  field  amplitude,  phase  and  gain: 

|a(t)|  *  £°.e(3'2)l'J^/2,  (11-76) 


<J)(X): 


f  -v/3 

i        -,  (11-77) 

\2)      2 


G(x)*i 

v  ;  9 


e(j/2)'/,V3T 


11-71 


The  gain  grows  exponentially  along  t  with  a  growth  rate 
proportional  to  j1/3. 

F.    BEAM  QUALITY 

Previously,  when  considering  the  interaction  of  the 
electron  and  optical  pulse  we  assumed  perfect  injection.  In 
reality  there  are  errors  involved.  The  two  factors  that  have 
the  most  influence  are  electrons  injected  off  axis  at  the 
initial  position,  y0,  and  electrons  injected  at  an  initial 

angle,  9y.  From  Reference  (2)  the  decrease  in  phase  velocity 

is 

AV  =  ^^(kVYo2  +  TV)-  (ii"79) 

1  +  K 
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The  change  in  phase  velocity  from  these  errors  is  made  more 
severe  with  more  undulator  periods  N.  A  longer  FEL  is  more 
sensitive  to  injection  errors. 

A  measure  of  the  quality  of  the  electron  beam  is  the 

emittance,  sy  =  y09y  ,  where  yc  is  the  rms  initial  position 

spread  of  electrons  along  y,  and  9y  is  the  rms  initial 
angular  spread  of  the  electrons  in  the  y  direction.  The  beam 
can  be   focused   to   minimize   either   y0   or   9y ,   however 

emittance  sy  remains  fixed.  The  minimum  change  in  phase 

velocity  occurs  when  the  beam  is  "matched"  so  that  both 
terms  in  (11-79)  contribute  equally.  Matching  occurs  when 
Kk0y0  «  y9y  .  For  a  matched  beam  Av  «  27iN0Kksy  /  y  . 

To  show  the  effects  of  beam  quality  we  need  to  look  at 
the  FEL  integral  equation  (11-73)  .  The  function  f  (q)  in 
equation  (11-73)  has  a  distribution  based  only  on  the 
initial  distribution  of  the  electrons  as  they  enter  the 
undulator.  Two  examples  of  distributions  are  those  from  the 
electron  phase  velocity  and  injection  angle.  The  spread  in 
phase  velocity  can  be  considered  a  Gaussian  distribution 

about  v0  with  a  standard  deviation  of  c?G  =  47iNAy  /  y  .   The 

angular  spread  leads  to  a  exponential  distribution  with  a 

standard   deviation   given   by   ae  =  4rtNy20y2  /(l  +  K2).   The 

distributions  are 
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fG(q)  = 

e"q  /2a< 

V2ucrG 

and 

■F    (rr\      — 

eq ' G" 

11-80) 


(11-81) 

The    characteristic    function   for   these    two    cases    are 

FG(q)   =   e-qV/2ac2  ,  (11-82) 

and 

Fe(q)   =   (11-83) 

1  -  iaex 

For  perfect  injection  aG=0  and  ae=0,  we  have  F(x)=l.  When  the 

beam   is  not  perfect   the   characteristic   function   |f(t)  I 

decays  with  a  characteristic  time  l/crG  or  l/ae.  F  (t)  depends 

on   the   shape   of   the   function   f  (q)  ,   and   different 
distributions  of  f (q)  effect  the  FEL ■ s  ability  to  bunch. 

G.     SHORT  PULSE  EFFECTS 

In  the  low  gain  FEL,  the  modes  of  light  in  the 
oscillator  can  develop  independently.  Since  the  electrons 
are  relativistic,  they  travel  to,  but  less  that,  the  speed 
of  light.  At  resonance,  one  wavelength  of  light  passes  over 
an  electron  as  the  electron  travels  through  one  undulator 
period.  Over  the  total  length  of  the  undulator  the  electron 

slides  behind  the  light  by  NA.,  called  the  slip  distance.  In 
the  transverse  dimension,  the  optical  mode  is  roughly  the 

27 


radius  of  the  optical  beam  which  is  (NA./ti)1/2.  These  modes 

evolve  as  determined  by  the  pendulum  equation  and  wave 
equation. 

In  order  to  follow  the  evolution  of  the  optical  field, 
we  must  set  up  multiple  sites  in  z  along  the  complex  wave 

envelope.  The  transformation  to  spatial  modes,  a— »a(z),  is 
the  equivalent  to  the  transformation  to  wave  numbers, 
a— »a(k)  .  The  longitudinal  position  on  the  z-axis  Z  is 
transformed  to  dimensionless  parameters  by  dividing  by  the 
slippage  distance  NX,  therefore  z— >Z/NX.  Througout  the 
remainder  of  this  paper  z  is  dimensionless  parameter. 

The   length   of   the   electron   pulse   le   becomes   the 

dimensionless  parameter  a2=le/NA..  The  pendulum  and  wave 
equations  of  motion  with  slippage  included  become 

L,=k|cOS(<;z_T+(|>z)  (11-84) 

and 

a*  =  -'h-^"x)  ■  (11-85) 

The  electron  current  is  jz  =  87re2p(z)NK2L2/y3mc2  ,  where  the 
electron  density,  p(z)  ,  is  the  density  at  the  site  z.  A 
site  in  the  electron  beam  interacts  with  and  amplifies  N 
optical  wavelengths,  while  the  optical  wave  passes  over  and 
attempts  to  bunch  electrons  in  a  slip  distance  along  the 
electron  beam. 
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Short  pulses  of  spontaneous  radiation  drive  a  typical 
FEL  oscillator.  When  the  length  of  the  electron  pulse,  le, 
is  comparable  to  the  slip  distance,  short  pulse  effects 
influence  the  interaction.  As  the  electron  and  optical 
pulses  pass  over  each  other,  the  pendulum  equation  (11-85) 
sees  a  continually  changing  optical  field,  az,  while  the 
wave  equation  (11-86)  sees  a  continually  changing  current 

density   j2.T   and   bunching   distribution  le'1^-')  .   In   the 


resonator,  the  optical  pulse  bounces  between  two  mirrors. 
The  distance  between  the  mirrors  is  S  which  is  necessarily 
greater  than  the  length  of  the  undulator,  L.  The  time 
interval  for  the  light-pulse  make  one  round  trip  is  2S/c. 

Desynchronism,  d,  is  the  difference  between  the  arrival 
times  of  the  electron  and  optical  pulses  at  the  beginning  of 

the  pass,  again  normalized  by  the  slip  distance  NA,.  If  d=0 
at  x=0,  the  electron  pulse  arrives  simultaneously  with  the 
light-pulse.  Moving  one  of  the  mirrors  slightly  by  AS 
changes  d  so  that  d=- AS/NX. 

H .    TAPER 

1.    Purpose  of  Taper 

In  a  successful  FEL  design,  the  optical  field  will 
eventually  saturate  due  to  over-bunching  of  the  electrons 
along  the  undulator,  thereby  changing  the  dynamics  of  the 
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electron  interaction.  As  the  electron  loses  energy  to  the 
optical  field,  the  phase  velocity  moves  off  resonance  by 

Av « 47tNA y/y « -2tc  .  One  way  to  extend  saturation  to  higher 
powers  is  to  taper  the  undulator  field.  Instead  of  a 
periodic  sinusoidal  magnetic  field  with  constant  amplitude, 
the  amplitude  of  the  field  is  either  increased  or  decreased 
progressively  along  the  -length  of  the  undulator.  Either 
shortening  the  undulator  period  along  the  undulator,  or 
changing  the  strength  of  the  magnetic  field  along  the 
undulator  tapers  the  undulator.  In  either  case,  acceleration 
of  the  phase  of  the  electron  results. 

This   acceleration   is   denoted  by   5,   or   the   taper 

parameter,  so  that  v  =  5,  and 

KAK 

5  =  -47tN ,  (11-86) 

1  +  K2 

where   the   undulator   field   is   tapered   by   AK  over   the 

undulator  length.  The  electron  pendulum  equation  becomes 

2;  =  v  =  5  +  |a|cos(<;  +  (j))  .  (11-87) 

While  the  optical  wave  equation  (11-86)  remains  unchanged. 

2 .    Taper  Phase  Space 

In  the  limit  of  low  gain  the  phase  space  paths  for  the 
tapered  undulator  are  given  by 

i 
v   =   [v02   -  26fc  -  C0)  -  2|a|[sinfc  +  <j>)  -  sinfe0   +  <f>)]]  ~2  (11-88) 
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The  stable  fixed  point  for  both  positive  (5>0)  and  negative 

taper  (5<0)  in  phase  space  is  at 

fou)=(cos"'(-5/|a|)-(t),0)  (11-89) 

The  unstable  fixed  points  depends  upon  the  direction  of  the 

taper.  The  unstable  fixed  point  for  6>0  in  phase-space  is 

located  at 

feu)  =  (27i-cos-1(-S/|a|)-(|),0)  (11-90) 

and  for  §<0  it  is  at 

feu)=(-cos-,(-S/|a|)-(|),0).  (11-91) 

The  formula  for  the  separatrix  is 

vs2  =25fcs-  Co  )+2|a|[sinfes +(}))- sinfe0  +<>)],  (11-92) 

where  ^^271-cos"1  ( -5/a)  for  5>0  and  ^0=-cos"1  ( -5/a)  for  8<0. 

Figure  II-3  shows  phase  space  plots  for  negative  and  postive 
tapers.  "U"  refers  to  the  unstable  fixed  point,  and  "S" 
refers  to  the  stable  fixed  point. 
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FEL  PHASE  SPACE  PATHS 
a=20       <|)=0         8=-2.57C 


-4Va 


FEL   PHASE    SPACE   PATHS 
a=20  (p=0  8=2 .  5JC 


I'.a 


-JC/2 


37C/2 


-JC/2 


371/2 


Figure  II-3  Phase  Space  Plots  for  Negative  and  Positive 
Taper 
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III.  TAPER  SEPARATRIX  HEIGHT 

To  find  the  height  of  the  separatrix  for  the  tapered 
undulator  it  is  convenient  to  make  a  few  variable  changes 
and  assume  that  there  is  no  initial  phase  (which  only  shifts 
the   separatrix   along   x)    C,   =  x,    v  =  y,    (j)  =  0,    and 

C,0  =  x0  .  Hence,  equation  (11-95)  becomes 

y2  =  28(x  -  x0)  +  2|a|(sin  x  -  sin  x0)  (III-l) 


-%/2 


37C/2 


Figure  III-l  Separatrix  for  negative  taper  with  q=0.4 

Figure  III-l  shows  the  separatrix  for  a  taper  8=-2.5  and  an 
optical  field  a0=20.  The  height  of  interest  is  the  maximum 
height  labeled  h  in  Figure  III-l.  To  find  the  maximum  value, 
we  need  to  take  the  derivative  of  equation   (III-l)  which 
gives 

2yy'  =  26  +  2|a|  cos  x  ,  (III -2) 
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where  y'  =  dy  /  dx  ,  and  set  the  slope  y'  equal  to  zero 


8  +  lal  cos  x 

y'  =  U =  0 

y 


III-3 


The  maximum  value  for  x  can  be  denoted  as  x  and  is  given  by 
8  +  |a|  cos  x  =  0    or    x  =  cos_1(-  q) ,  (III-4) 

where  q  =  8  /|a|  .  Note  that  q  is  restricted  to  values  between 
1  and  -1. 

A.     POSITIVE  TAPER 

For  S>0,  x0  =  2tt-cos"1  ( -q)  and  substituting  x  for  x  gives 


y2  =  2S[cos_1(-  q)  -  (2tt  -  cos_1(-  q))]  + 

2|a|[sin(cos_1(-  q)j  -  sin(27i  -  cos_1(-  q))j 


(III-5) 


Using  the  identity  sin(27i-x)  = -sin(x)  we  can  simplify  equation 
(III-5)  to  get 

y2    =   28[2cos_1(-  q)  -  2tt]  +  2|a|[2  sin(cos_1(-  q))]  .  (III-6) 

Factoring   out    |a|    we    have 

y2    =   4|a|[q(cos_1(-  q)  -  n)  +  sin(cos_1(-  q))]  .  (Ill -7) 

Using     the     identities     cos"'(-  x)  -  n   =   -  cos_1(x)     and 

sin(cos~^(-  q)J  =  -yjl  —   q2  ,  we  can  consolidate  equation  (III-7) 
to  get 

y2  =  4|a|-  q  cos":(q)  +  Jl  -   q2   .  (III-8) 
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Taking  the  square  root  of  both  sides,  we  get  the  maximum 
value  of  y  to  be 


y 


=  \2^M\~   ^  cos-i(q)  +  V1  -  q\ 


,  i 

2 
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The  height  of  the  separatrix,  h  is  twice  the  value  of  y 
therefore , 


h  =  4Ja 


■,  i 


-  q  cos_1(q)  +  -y/l  ~  q' 


(III -10) 


B.    NEGATIVE  TAPER 

For  negative   q,   x0  =  -  cos_i(-  q)  .   Therefore   equation 
(III-l)  becomes 

y2  =  2§[cos":(-  q)  -  (-  cos_1(-  q))l  + 

\        I  \  t  M  (III-ll) 

2|a|[sin(cos"1(-  q)j  -  sin(-  cos  1(-   q))J  . 

Utilizing   the    identity    sin(-  x)  =  -sin(x),    it   becomes 

y2    =   4§[cos_1(-  q)]  +  2|a|[-  sin(cos_1(-  q))  -  sin(cos_1(-  q))]       (111-12) 


Consolidating  and  factoring  out  |a|  gives 

y2  =  4|a|[q(cos_1(-  q))  +  sin(cos_1(-  q))j  . 

Take  the  square  root  of  both  sides  gives 


y  =  2^J\a\q  cos  1(-  q)  +  V1  -  ^ 


and   then   the   height    is 


h   =  4Ja 


q 


cos'1^-  q)  +  V1  ~ 


i 

2 


(111-13) 


111-14) 


111-15 
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Equations  (111-10)  and  (111-15)  are  equivalent  with  the  |q 

Therefore  we  can  say  for  all  values  of  q 

,  i 


h  =  4,/ a 


-  |q|  cos_1(jq|)  +  V1  -  q" 


(111-16) 


C.     POSTIVE  AND  NEGATIVE  TAPER 

Since  the  separtrix  height  for  the  non-tapered  case 
(5=0)  is 


h  =  4Ja 


|q|  cos  \\q 


l)  +  VT^ 


111-17' 


Figure  III-2  shows  a  graph  of  equation  (111-18)  . 


S/R 


0.5- 


0  .  5 
o/a 


Figure  III -2  Change  in  Separatrix  Height  8/a 
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D.    LINEAR  APPROXIMATION 

A  Taylor  expansion  inside  the  brackets  of  equation 
(111-18)  and  using  the  expansion  of  arccosine 


71 


cos  1{q)   = q  +  .  . 

2 


111-18) 


gives 


h  =  4Ja 


71 

a7-q 


1   2 
+  l--q- 


(111-19 


which  simplifies  to 


h  =  4Ja 


7T       1   • 

1  —  q  +  -  q 

2     2 


111-20 


h  =  4Ja 


t  71 

1  -  -  q 


V 


( 


=   4,  a 


4    ; 


A 


1  -  0.785  - 
ay 


111-21) 


From  equation  (111-21)  it  can  be  shown  that,  as  8  — >  0  the 
separatrix  height  h  — >  4^|a|  ,  or  the  size  becomes  the  same  as 
the  untapered  case.  As  the  optical  field  a  — >  oo  and  the 
height  h  — >  4j|a|  the  same  as  5=0.  As  the  optical  field 
becomes  zero  a  — >   0  then  the  height  h  — >  0 . 
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IV.   PARAMETERS  FOR  USE  IN  SIMPLE  MODEL  OF  A  HIGH  GAIN  FEL 

To  develop  a  simple  model  for  a  high  gain  FEL  we  need 
to  first  take  a  close  look  at  the  parameters  that  dominate 
the  electron  interaction  with  the  optical  field  in  the  high 
gain  regime.  Table  IV- 1  is  a  table  of  several  proposed  and 
well -documented  high  gain  FELs . 


Units 

LCLS 

TESLA 

ELF85 

ELF86 

ELF87 

MIT 

KE 

MeV 

1500 

35000 

3.3 

3.5 

3.5 

2.3 

I 

amps 

3400 

5000 

500 

850 

1000 

930 

X 

cm 

1.5xl0"8 

lxlO"8 

0.76 

0.72 

0.175 

0.125 

L 

m 

100 

60 

3 

3 

4 

1.5 

J 

- 

27200 

35000 

3300 

5000 

5600 

3800 

°x 

- 

0.45 

0.003 

0.045 

0.04 

0.12 

0.16 

a0 

- 

0 

0 

12 

14 

0.34 

0 

^o 

cm 

3 

5 

9.8 

9.8 

9.8 

3.14 

N 

- 

3300 

1200 

30 

30 

40 

48 

Np 

- 

0.3 

0.4 

11 

11 

5.7 

2.3 

<*G 

- 

8.4 

3.7 

12 

12 

32 

3.6 

1. 

cm 

0.008 

20 

450 

1498 

300 

900 

NX 

cm 

3xl0~5 

2xl0"5 

23 

22 

7 

6 

Slip 

- 

186 

1000 

19 

68 

42 

6,000 

W0x 

cm 

0.002 

0.006 

5 

5 

1.6 

0.8 

W0y 

cm 

0.002 

0.006 

1.5 

1.5 

1.5 

0.8 

re 

cm 

0.003 

.0018 

0.4 

0.3 

0.6/    1.2 

0.4 

s 

cm-rad 

0.00015 

0.0003 

0.37 

0.07 

0.07 

0.16 

F 

- 

1 

1 

0.026 

0.024 

0.24 

0.28 

Table  IV- 1  Table  of  High  Gain  FEL  Parameters 


39 


The  following  chapter  is  a  description  and  analysis  of 
each  FEL  in  order  to  characterize  them  with  our  simple 
model . 

A.     DESCRIPTION  OF  LCLS  AT  SLAC 

The  Linac  Coherent  Light  Source  (LCLS)  proposal  is  to 
use  the  last  kilometer  of  Stanford  Linear  Accelerator  Center 
(SLAC)  to  build  a  FEL  able  to  achieve  a  1.5  A  x-ray 
wavelength.  An  electron  gun  will  inject  an  electron  bunch 
that  will  be  accelerated  along  the  last  kilometer  of  the 
SLAC  accelerator  to  14.35  GeV  energy  with  a  peak  current  of 

3400  A.  The  normalized  emittance  is  eN  of  1.5  mm-rad.  After 

the  undulator,  the  electron  beam  will  be  dumped.  The 
accelerator   is   tunable   from   14.3   to   4.5   GeV   which 

corresponds  to  x-ray  wavelengths  of  1.5-15  A.  The  electron 
beam  particle  density  is  p=2.5xiols  cm"3  in  a  beam  radius  of 

re=30u.  The  undulator  will  be  a  linearly  polarized  with  a 

length  of  L=112  m,  period  A,0=3 . 0  cm,  number  of  periods 
N=3328,  peak  magnetic  field  strength  B=1.32T  for  a  resultant 
undulator  parameter  of  K=2.6,  and  a  dimensionless  current  of 
j=27, 000 . 

The  FEL  will  be  a  single-pass  amplifier  starting  from 
noise  or  spontaneous  emission;  a  process  called  Stimulated 
Amplification  from  Spontaneous  Emission  (SASE) .  Due  to  the 
lack  of  mirrors  at  this  wavelength  , numerous  proposals  are 
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being  investigated  to  guide  the  light  once  it  leaves  the 
FEL. 

B.  TESLA  FEL  AT  DESY 

The  Deutsches  Elektronen-Synchrotron  (DESY)  in  Hamburg, 
Germany,  is  also  planning  to  develop  an  x-ray  FEL  using  the 
TeV  Energy  Superconducting  Linear  Accelerator  (TESLA) .  The 
electron  beam  will  be  supplied  by  two  16  km  long 
superconducting  linear  accelerators  which  were  initially 
designed  to  collide  electrons  and  positrons. 

The  specifications  for  the  TELSA  FEL  are  still  under 
consideration  and  the  project  is  being  planned  in 
conjunction  with  the  LCLS  with  a  later  timeline.  Many 
parameters  will  be  based  on  experiments  conducted  in 
preparation  for  LCLS.  The  projections  are  based  on  beam 
energy  of  25GeV  and  a  peak  current  of  5kA. 

C.  ANALYSIS  OF  SLAC  AND  TESLA  PARAMETERS 
1 .    Wavelength 

The  purpose  of  these  two  lasers  is  to  create  a  source 
of  high  power,  coherent,  hard  x-rays  at  a  1  A  wavelength  for 
numerous  experimental  uses.  The  lack  of  normal  incidence 
optics  at  this  wavelength  creates  many  unique  problems.  Only 
grazing-incidence  optics  are  possible,  and  they  are  still 
under  development.  The  limited  optics  precludes  the  use  of 
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an  oscillator  configuration  and  requires  the  FEL  to  operate 
as  a  single  pass  amplifier. 

2 .    Reducing  j 

The  dimensionless  current  j  for  both  x-ray  lasers  is  a 
large  number.  There  are  several  factors,  which  lower  the 
effective  j  for  actual  use  in  our  simulations. 


a)         Filling  Factor 

In  a  typical  FEL,  the  radius  of  the  electron  beam 
is  smaller  than  that  the  optical  mode.  The  filling  factor  is 
the  cross  sectional  area  of  the  electron  beam  divided  by  the 
cross  sectional  area  of  the  optical  beam.  For  LCLS  and  TESLA 
the  cross  sectional  area  of  the  electron  beam  is  comparable 

to  the  cross  sectional  area  of  the  optical  beam,  so  F«l .  For 

LCLS,  the  radius  of  the  electron  beam  re  is  about  0.03  mm= 

3  0  urn  so  the  normalized  pulse  width  ax~0.45.  The  radius  of 

the  optical  field  is  approximately  w0«0.02  m=  2  urn. 


Gaussian 
OpicalBeam 


Figure  IV- 1  Gaussian  X-ray  Beam  and  Top-Hat  Electron  Beam 
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b)         Betatron   Oscillation 

The  number  of  electron  betatron  oscillations  along 
the  undulator  in  the  x-ray  FELs  is  small.  The  number  of 
betatron  oscillations  along  the  undulator  is  inversely- 
proportional  to  the  energy  of  the  electron  beam.  If  there 
were  a  large  number  of  betatron  oscillations,  then  the 
average  amount  of  time  electron  beam  spend  within,  the. 
optical  field  is  of  little  concern.  Betatron  oscillations 
cause  the  electrons  to  pass  from  the  outside  of  the  optical 
mode,  where  the  fields  are  weak  to  the  center  of  the  mode, 
where  the  fields  are  strong.  In  order  to  account  for  the 
reduced  bunching  drive  of  the  electrons  due  to  betatron 
oscillations,  due  to  the  dimensionless  current  j  is 
decreased  by  half. 

3 .    Rayleigh  Length 

The  waist  size  of  the  x-ray  beam  is  very  small  with  a 

w0=20   urn  radius.   A  small  Rayleigh  waist  means  a  small 

Rayleigh  length  and  large  diffraction  pattern  a  large 
distance  away.  The  x-ray  beam  diffracts  away  from  the 
electron  beam  contributing  to  reduced  coupling.  The  Rayleigh 

length,   z0  =  7rw0~  /  k  ,   for  the  LCLS  this  is  approximately 

30m,  about  one-third  the  length  of  the  undulator.  In  a  high- 
gain  FEL  such  as  LCLS  the  gain  is  exponential  along  the 
undulator.  For  LCLS,  the  gain  in  a  Rayleigh  length  is  about 
100m  as  shown  in  Figure  IV-2. 
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electron 
beam 


optical  beam  spreading 


Figure  IV- 2  Discrete  comparison  of  gain  versus  diffraction 

To  illustrate  the  reduced  role  of  diffraction  of  the  x- 
ray  beam,  divide  the  undulator  into  three  discrete  pieces. 
At  the  beginning  of  the  undulator  the  optical  power  is  taken 
to  be  1  as  shown  in  Figure  IV-2  on  the  left.  In  one  Raleigh 
length,  z0,  along  the  length  of  the  undulator,  the  gain  is 
100,  but  only  within  the  region  of  the  electron  beam.  The 
power  in  the  center  is  100,  while  the  x-ray  beam  outside  the 
electron  beam  has  had  no  gain  and  remains  power  1 .  The 
strength  of  the  diffracted  x-ray  beam  is  only  1%  of  the 
power  of  the  x-ray  beam  in  the  center.  In  the  second 
Rayleigh  length,  the  x-ray  within  the  electron  beam  region 
has  been  amplified  another  factor  of  100.  The  power  in  the 
center  x-ray  beam  after  two  Raleigh  lengths  is  1002.  The 
power  of  the  diffracted  beam  was  amplified  to  100  in  the 
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first  Rayleigh  length  and  remains  100.  Again  this  only  1%  of 
amplified  x-ray  beam  within  the  electron  beam  region.  In  the 
third  length,  the  power  at  the  end  of  the  undulator  has 
become  1003,  while  the  diffracted  light  has  only  power  of 
1002.  The  power  of  the  diffracted  region  is  again  only  1%  of 
the  power  of  the  beam  in  the  center  at  the  end  of  the 
undulator. 

In  reality  this  is  a  continuous  process  along  the 
undulator.  But,  just  as  for  this  discrete  example,  the 
continuous  gain  within  the  central  area  of  the  electron  beam 
far  out -weighs  the  gain  at  the  edges  of  the  electron  beam. 
To  account  for  this  reduction  in  coupling,  the  value  of  j 
can  be  reduced  by  half,  so  that  considering  betatron  motion 

and   diffraction   the   final   value   is   j«6700.     This   is 

comparable  to  the  other  long  wavelength  FELs  in  Table  IV- 1. 

4 .    Emittance 

The  emittance  for  the  high  gain  FEL  can  be  extremely 
important.  The  x-ray  field  growth  is  exponential  with  an  e- 

folding  time  of  x^  »  2/(j/2)1/3V3  »  0.15  for  j=5000.  For  LSCS 
the  emittance  of  Ge=0.7,  leads  to  a  e-folding  decay  time  of 
xa  ~  1  /  ae  «  1.4  .  Since  the  decay  time  caused  by  emittance 
is  greater  than  the  growth  e-folding  time  driven  by  the 
electron  beam,  ta  «  1.4  >  t^  «  0.1.  The  net  growth  rate  of 
the  x-ray  beam  is  significantly  reduced  along  the  undulator. 


45 


5.  Slip 

The  electron  pulse  length  is  about  186  time  longer  than 

the  slip  distance  NA,.  The  extremely  short  x-ray  wavelength 

makes  the  slippage  distance  very  small  and  short  pulse 
effects  are  non-existent.  Each  part  of  the  x-ray  pulse 
experiences  gain,  which  is  proportional  to  the  local  density 
of  electrons. 

6.  Initial  Optical  Field 

Since  no  other  laser  can  radiate  near  x-ray 
wavelengths,  the  x-ray  FEL  must  seed  itself.  The  initial  x- 
ray  field  a0  for  an  x-ray  laser  is  necessarily  zero,  and  the 
laser  depends  on  the  SASE  (stimulated  amplification  of 
spontaneous  emission)  to  begin.  A  simple  simulation  model  is 
to  start  with  a  very  small  initial  optical  field  that  is 
representative  of  the  spontaneous  emission  created  in  the 
beginning  section  of  the  undulator.  In  the  first  portion  of 
the  undulator,  the  model  assumes  that  an  initial  optical 
field  starts  from  random  fluctuations  in  the  electron  beam 
density.  The  x-ray  mode  in  SASE  is  known  to  be  linear, 
starting  from  random  noise,  and  growing  through  mode 
competition  to  the  desired  wavelength.  Since  the  field 
growth  in  the  beginning  of  the  undulator  is  linear,  instead 
of  the  exponential  gain  seen  later  along  the  undulator,  a 
slightly  longer  undulator  is  required. 


46 


D.     ELECTRON  LASER  FACILITY  (ELF) 

The  Electron  Laser  Facility  (ELF)  is  an  offshoot  of  the 
Experimental  Test  Accelerator  (ETA)  located  at  Lawrence 
Livermore  National  Laboratory.  The  ELF  was  built  to  utilize 
the  electron  beam  accelerated  by  ETA.  In  each  year,  1985, 
1986  and  1987  a  paper  was  published  updating  that  years 
experiments  at  ELF  (references  4,5,  and  6)  .  The  main 
components  of  the  laser  remain  essentially  the  same  with 
several  modifications  from  year  to  year.  Here,  the  three 
phases  of  the  ELF  are  referred  to  as  ELF85,  ELF86,  ELF87. 
ELF8  5  had  beam  energy  of  3.3  MeV  and  ELF8  6  and  ELF8  7  had  a 
beam  energy  of  3.5  MeV.  ELF85  and  ELF86  used  a  3m  long 
undulator  which  was  extended  in  ELF87  to  4  meters.  The 
undulator  for  all  three  is  composed  of  specially  shaped 
solenoids  with  a  peak  magnetic  field  of  5  kG.  The  wiggler 

period  remained  A,0=9.8  cm  for  all  three.   The  number  of 

periods  was  N=30  and  N=40  respectively.  The  peak  current 
initially  was  500  A  for  ELF85,  then  increased  to  850  A  for 
ELF86,  and  finally  to  1000  A  for  ELF87.    The  wavelengths 

produced  were  A.=  0.76  cm,   0.72  cm,  and  0.175  cm.  For  each 

experiment,  the  optical  field  was  contained  in  a  2.9  cm  by 
9 . 8  cm  rectangular  waveguide  and  the  FELs  were  operated  as 
single  pass  amplifiers. 
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E. 


ELF85 


ELF85  has  a  dimensionless  current  j=33  00.  The 
dimensionless  beam  size  ax  is  0.045,  and  the  initial  optical 
field  of  a0=12  was  produce  by  a  60kW  pulsed  magnetron. 

1.    Electron  Beam  to  Optical  Field  Ratio  the  Filling 
Factor 

A  sketch  of  the  cross  section  of  electron  beam  and 
optical  mode  is  shown  in  Figure  IV- 3  and  the  filling  factor 
is  about  F=0.03. 
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Figure  IV-3  ELF  85  optical  field (white) ,  electron  beam(gray) 

diagram 

The  electron  beam  is  small  enough  so  that  electrons  see  only 
the  peak  of  the  optical  beam.  There  are  a  large  number  of 
betatron  oscillations  along  the  undulator,  unlike  LCLS,  but 
the  electrons  continue  to  see  the  peak  optical  field.  The 
betatron  oscillations  increase  little  effect,  since  they 
remain  within  the  optical  field.  There  is  some  radiation  in 
the  harmonics  due  to  the  oscillation,  but  not  enough  to  be 
significant . 
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2.  Slip 

The  electron  pulse  is  20  slippage  distances  in  length 
for  ELF85.  This  is  much  smaller  than  for  the  LCLS ,  but  still 
much  larger  than  1.  Since  the  slip  is  larger  than  1,  short - 
pulse  effects  can   still  be  neglected. 

3 .  Emittance 

The  emittance  for  the  beam  was  quite  large  e=0.37  cm- 
rad.  This  large  emittance,  assuming  a  matched  beam,  means  a 
ae=22,   and   therefore   a   characteristic   decay   time   of 

Ta=l/ae=0 .  045 .  Gain  degradation  due  to  the  emittance  is  quite 
significant . 

F.    ELF  86 

The  major  difference  between  the  ELF85  and  ELF86  is  the 
modification  in  the  Experimental  Test  Accelerator  (ETA)  to 
use  a  field  emission  cathode.  This  increases  the  peak 
current  and  makes  the  pulse  longer  and  the  emittance 
smaller.  The  undulator  dimensions  remained  the  same. 
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1.    Electron  Beam  to  Optical  Field  Ratio,  the  Filling 
Factor 


Figure  IV-4  ELF  86  optical  field,  electron  beam  diagram 

The  ratio  of  the  electron  beam  area  to  the  optical  beam 
area  decreased  slightly  giving  a  filling  factor  F=0.024,  as 
seen  in  Figure  IV-4.  The  electron  pulse  length  remained  at 
20  slippage  distances  with  minimal  short  pulse  effects.  The 

dimensionless  beam  size  ax  remained  nearly  the  same  due  to  a 

slightly  smaller  optical  wavelength  0.72  cm.  The 
dimensionless  current  increased  to  j=5000.  The  input  optical 
field  is  provided  by  a  lOOkW,  35GHz  magnetron  giving  a0=14 . 

The  emittance  was  significantly  reduced  to  s=0.07  cm- 

rad.  The  overall  increase  in  current  and  better  beam  quality 
has  increase  the  power  output  from  80  MW  for  ELF85  to  180MW 

for  ELF86  and  the  efficiency  increased  from  n  =  5%  to  r|=6%. 

G.     ELF87 

A  major  change  in  ELF87  was  the  use  of  a  high-current 
induction  linac.  This  linac  provided  a  beam  energy  of  3.5 
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MeV  and  a  peak  current  of  1  kA.  The  dimensionless  current  j= 
10  0  0  and  the  undulator  was  also  extended  to  L=4  m. 

Another  major  change  was  the  shape  of  the  electron 
pulse.  The  beam  length  cross  section  was  oval  as  seen  in 
Figure  IV-5.  However,  this  beam  shape  has  little  effect  on 
the  model. 

The  shift  in  the  optical  wavelength  to  ^=0.175.  cm 
causes  the  optical  mode  to  become  much  smaller.  The  ratio  of 
electron  beam  area  to  optical  mode  area  is  much  smaller.  A 
process  called  optical  guiding  also  causes  the  optical  mode 
to  remain  small. 
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Figure  IV-5ELF  87  optical  field,  electron  beam  diagram 

The  net  effect  is  a  larger  filling  factor,  F=0.24,  and  more 
interaction  between  the  light  and  the  electron  pulse. 


H. 


MIT  FEL 


The  MIT  FEL  is  operated  by  the  Department  of  Physics 
and   Research   Laboratory   of   Electronics,   Massachusetts 
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Institute  of  Technology.  This  FEL  was  designed  to  operate  at 
millimeter  and  sub-millimeter  wavelengths.  The  MIT  FEL  has 
not  been  studied  as  closely,  since  it  has  lower  gain  than 
ELF.  The  most  significant  difference  between  this  laser  and 
the  ELF  laser  is  the  undulator.  The  MIT  laser  uses  a  helical 
undulator . 

1.  Natural  Focusing 

An  advantage  of  the  helical  undulator  is  the  natural 
focusing  of  the  electron  beam.  The  magnetic  field  is  always 
pushing  the  electrons  toward  the  center  of  the  beam.  Far  an 
x-ray  FEL,  any  added  benefit  from  natural  focusing  could 
result  in  a  significant  increase  in  gain. 

2 .  SASE 

The  other  significant  feature  of  the  MIT  laser  is  the 
SASE  startup  at  long  wavelengths.  From  its  operation,  it  can 
be  seen  that  SASE  occurs  within  about  the  first  20  undulator 
periods.  After  the  first  20  periods,  the  FEL  field  growth 
becomes  exponential  and  can  be  simulated  using  a  small 
initial  optical  field  a0 . 

3 .  Electron  Beam 

A  multielectron  gun  with  a  Pulserad  110A  accelerator 
creates  the  electron  beam.  The  beam  is  of  relatively  low 
kinetic  energy,  2.3  MeV,  with  a  peak  current  of  930  amps. 
The  pulse  length  is  7200  times  the  slippage  distance,   so 
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that  the  light  pulse  interacts  with  the  electrons  along  the 
entire  length  of  the  undulator. 

I .    CONCLUSION 

In  careful  study  of  these  FEL's,  one  can  discern  which 
parameters  have  the  most  influence  on  the  growth  rate  and 
FEL  interaction.  A  simulation  model  can  often  be  simplified 
by  concentrating  on  those  parameters  which  create  the  large 
growth  rate  or  decays,  and  ignoring  those  which  have  no 
discernible  effect  on  the  electron/optical  field 
interaction.  The  parameters  developed  in  this  chapter  can  be 
used  in  a  simple  model.  This  is  left  for  further  study. 
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V.    DARMSTADT  FEL 

The  next  section  is  an  analysis  of  the  Darmstadt  FEL 
utilizing   computer   simulations.   The   Darmstadt   FEL   is   a 

medium  gain  laser  and  operates  between  X=6.G     and  7 .  8um 
wavelength,  with  electron  beam  energies  between  32  and  29 
MeV.   It   is   driven   by   the   S-DALINAC,   a   re-circulating 
electron  linac  using  superconducting  acceleration  cavities, 
which  allow  for  a  continuous  wave  electron  beam. 


A. 


DARMSTADT  PARAMETERS 


The  parameters  for  the  Darmstadt  FEL  are  shown  in  Table 
V-l.  The  primary  use  for  the  Darmstadt  FEL  is  to  provide 
infrared  light  for  medical  research.  Most  of-  the 
experiments  involve  medical  uses  such  as  ablation  of 
biological  tissue. 


Electron  Energy-KE 

Peak  Current -I 

Electron  Pulse  Length-le 

Number  of  Periods -N 

Undulator  period-A,0 
Undulator  Parameter-K 

Dimensionless  current -j 

Energy  Spread-AE/E 
Q-f actor 


31  MeV 
2.7  A 

6  0  0um 

80 

32mm 
0.78 

1.6 
0.0023 

110 


Table  V-l  Darmstadt  parameters 
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B.     BEAM  SIZE  EFFECTS  ON  GAIN 

Initially  the  parameter  values  obtained  for  the 
Darmstadt  FEL  were  not  properly  understood.  The  electron 
beam  radius  was  given  at  the  end  of  the  undulator,  and  it 
was  mistaken  to  be  the  radius  at  the  center  of  the 
undulator.  After  several  communications  and  some  analysis, 
the  discrepancy  was  discovered.  As  a  result  of  angular 
spreading  along  the  undulator  the  beam  radius  was  0 . 2mm  at 
the  waist,  and  a  little  less  than  1mm  at  the  end  of  the 
undulator.  For  an  electron  beam  with  Darmstadt's  parameters 
and  an  electron  beam  waist   radius   re=lmm,   the   angular 

distribution  would  be  gq=7.3.       In   the   time   it   took   to 

discover  the  error,  an  analysis  of  the  influence  of  angular 
spread  on  gain  was  conducted.  Despite  not  having  any 
relevance  to  the  Darmstadt  FEL,  the  analysis  was  interesting 
and  may  have  future  relevance.  Therefore  the  results  are 
included  here.  The  analysis  below  was  conducted  assuming  an 
electron  beam  waist  radius  of  1mm. 
1 .    Theoretical  Background 

An  angular  distribution  of  ae=7.3  is  significant  enough 

to  degrade  the  gain  of  an  FEL.  A  large  electron  beam  radius 
means  that  there  are  a  significant  number  of  the  electrons 
which  are  injected  into  the  undulator  that  are  off -set  from 
the  axis  of  propagation  and  at  an  angle  with  respect  to  that 
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axis.  The  larger  the  diameter,  the  more  electrons  that  are 
injected  off  center  and  at  larger  angles.  The  electrons  on 
the  outer  portion  of  the  beam  experience  a  larger  magnetic 
field.  These  outer  electrons  are  "wiggled"  harder  and 
therefore  travel  on  a  slightly  longer  path.  As  a  result  of 
the  longer  path,  the  electron  slows  down  in  the  longitudinal 
direction  with  respect  to  an  electron  in  the  center  of  the 
beam.  The  spread  in  longitudinal  velocities  creates  a  spread 
in  phase  velocity.  A  spread  in  phase  velocity  then  results 
in  a  reduction  in  gain. 

The  same  argument  can  be  used  for  electrons  that  are 
injected  at  an  angle  with  the  axis.  An  electron  injected  at 
an  angle  moves  toward  (or  away  from)  the  more  intense 
magnetic  field  off  axis.  The  electron  is  pushed  back  to  the 
center  harder  than  the  other  electrons  and  travels  a  longer 
distance.  The  longer  distance  means  a  greater  spread  in 
phase  velocity  and  less  gain. 

2 .    Method  of  Analysis 

In  order  to  determine  the  effect  of  beam  size  on  the 
gain,  the  effect  of  value  of  gq  on  the  gain  is  analyzed.  The 
gain   for  each  ae  must   be   calculated  at   the   value   of 

desynchronism  giving  maximum  gain.  For  each  ae,  there  is  a 
different  desynchronism  value  d  for  which  gain  is  a  maximum. 
The  amount  of  time  to  run  enough  simulations  to  determine 

the  maximum  gains  for  each  gq   would  have  been  prohibitive. 
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Gain  versus  desynchronism  curves  were  plotted  for  ae=7.3  and 
CTe=3 ,   as   shown   in   Figure   V-l.   In   order   to   use   this 

information  to  determine  how  the  maximum  d  depends  on  ae, 
linear  interpolation  was  used.  Assume  a  linear  relationship 
between  the  maximum  desynchronism  d  and  ae  using  the  form 
d  =  ma6  +  b  for  a  straight  line.  The  maximum  gain  for  each 
value  of  ae  can  be  estimated  in  subsequent  simulations. 

Maximum  gain  was  found  for  Ge=7.3  at  d=0.003,  and  for  Ge=3 
at  d=0.0025.  After  a  small  amount  of  algebra  we  find  that 
m=-0. 00058   and  b=0.0066.   To  test   the  validity  of   this 

estimation,  Ge=5  was  estimated  to  be  a  maximum  at  d=0.0037. 

A  series  of  simulations  shows  that  in  fact,  the  maximum 
occurred  at  0.0036;  in  good  agreement  with  the  interpolated 
formula.  In  Figure  V-l  it  can  be  seen  that  there  is  a  flat 
plateau  about  the  maximum  of  the  curve,  so  our  approximation 
for  the  maximum  is  accurate  enough.  In 

Figure  V-2  the  gain  found  for  each  ae,  at  the  maximum  d 

determined  by  the  interpolation  formula  above,  is  plotted.  A 

value  of  ae=0  means  perfect   injection,   no  deviation  in 

either  angle  or  position.   The  actual  beam  radius  re=0.2mm 

corresponds  to  ae=0.3  which  is  nearly  perfect.  The  gain  for 

ae=0.5  is  about  3%  less  than  for  perfect  injection  which  is 
an  insignificant  reduction. 
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Figure  V-l  Gain  Versus  De synchronism 

On  the  other  hand,  the  original  value  of  Ge=7.3  results  in  a 

77%  decrease  in  gain  which  is  not  negligible.  In  fact,  it 
would  mean  the  Darmstadt  FEL  would  not  work. 

From  Figure  V-2  the  effects  of  an  angular  or  radial 
spread  on  the  reduction  of  Darmstadt  gain  can  be  straight- 
forwardly determined. 
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Figure  V-2  Oq   versus  Gain 

The  final  power  P  versus  desynchronism  d  curves  for  <Jq=1  .3 

and  <7e=3  were  plotted,  as  seen  in  Figure  V-3 .  The  final 
power  in  strong  fields  at  saturation  peaks  at  small  values 
of  d  where  gain  is  small,  and  extends  out  to  larger  values 
of  d. 
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Figure  V-3  Power  Versus  De synchronism 

C.    TAPER  EFFECTS 

Using  the  actual  Darmstadt  FEL  electron  beam  radius  of 
re=0.2mm,  and  the  parameters  from  Table  V-l,  the  power 
versus  desynchronism  graphs  for  several  tapers  were 
determined.  The  following  positive  and  negative  tapers  were 

analyzed:  8=-5.47r,  -2.5ti,  -0.6371,  0,  1.3ti,  2.57t,  and  5.471. 

1.    Compiled  Data 

Figure  V-4  and  Figure  V-5  are  3-D  graphs  that  depict 
power  (on  the  vertical  axis)  as  a  function  of  both  taper  6 
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and  desynchronism  d.  The  two  figures  contain  the  same  data 
but  presented  iron  different  points  of  view. 
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Figure  V-4  3-D  Plot  of  Power  Versus,  d,  and  6 (Positive  View) 

The  power  versus  d  curves  ?(d)  for  each  value  of  taper 
5  have  similar  structures.  When  d=0,  there  is  no  power 
output  at  any  6.  Ac  approximately  d=0.008,  the  output  power 

is  maximum  for  each  taper  6.  The  de synchronisms  in  the  range 
0<d<C.0C8  dc  not  produce  any  power  for  any  taper. 
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Figure  V-5  3-D  Plot  of  Power,  d,  and  6  (Negative  View) 

The  sharp  rise  in  power  and  high  peak  power  is  typical  of 

many  FELs  operating  in  "his  parameter  regime.  The  largest 

peak  power  occurs  at  5=0. 

A  high  quality  factor  Q  means  chere  are  small  optical 
loses  per  pass  through  the  undulator.  These  losses  include 
mirror  loss  and  mirror  transparency  where  optical  power  is 
removed  for  the  end  user.  A  high  Q  leads  to  sharper  features 
in  the  desynchronism  curves 
explained  in  the  next  secti< 
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2 .    Untapered  Power  Curve  Analysis 

To  show   more  detail  of  the  data  in  Figure  V-4  and 
Figure  V-5,   Figure  V-6  and  Figure  V-9  show  power  versus 

desynchronism  with  no  taper  8=0  and  the  highest  value  of 
negative  taper  8=-5.47C. 
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Figure  V-6  P  Versus  Desynchronism  6=0 
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In  Figure  V-6,  the  FEL  exhibits  limit -cycle  behavior 
for  desynchronisms  between  d=0.0008  and  d=0.011  (ref  (2)), 
and  thus  is  not  stable.  Since  the  power  levels  are  not 
stable  the  power  in  this  region  is  an  average  of  the  power 
oscillations . 

3 .    Simulation 

Figure  V-7  is  the  output  from  a  single  simulation  run 
for  the  Darmstadt  FEL.  The  input  variables  for  this  run  are 
displayed  in  the  top  square.  There  are  nine  output  windows 
displayed  below  the  input  window. 

The  top-left  window  shows  the  dimensionless  optical 

field  amplitude   |a(z)  I  at  the  end  of  n=6000  passes.  The 

middle -left   window   shows   the   evolution   of   the   field 

amplitude  after  at  each  pass   |a(z,n)  I,   in  this  case  for 

n=6000  passes.  The  lighter  the  color,  the  higher  the 
amplitude,  as  is  the  case  for  each  middle  window.  In  the 

bottom-left  window,  the  parabolic  electron  pulse  j  (z-x)  is 

plotted  for  reference  at  the  beginning  and  end  of  the 
undulator;  their  separation  shows  the  slippage  along  the 
undulator . 

The  bottom-center  window  plots  for  reference  the  weak- 
field  gain  spectrum  as  a  function  of  phase  velocity,  v.  The 
middle-center  window  shows  the  optical  power  evolution 
P(v,n)  as  a  function  of  phase  velocity  v  at  each  pass  n.  The 
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top-center  window  show  the  final  optical  power  spectrum  P(v) 
at  the  end  of  n=6000  passes. 

The  bottom-right  window  shows  the  total  optical  power 
evolution  ?  (n)  as  a  function  of  the  pass  number  n.  The 
center-right  window  shows  the  evolution  of  the  electron 

distribution  in  phase  velocity  f(v,n)  as  a  function  of 
phase  velocity  and  number  of  passes  n.  The  upper- right 
window   displays   the   final   electron   phase   velocity 

distribution  f  (v,6000)  at  the  end  of  n=600  0  passes. 
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Figure  V-7    Trapped   Particle    Instability 
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4.    Limit-Cycle  Behavior  in  the  Untapered  case 

In  the  untapered  undulator,  the  trapped-particle 
instability  occurs  for  small  values  of  desynchronism  at 
d=0.012.  In  strong  optical  fields,  electrons  become  trapped 
in  the  deep  potential  wells  in  phase  space.  These  trapped 
electrons  oscillate  and  radiate  into  side-  bands  off  the 
main  wavelength  of  the  optical  field.  In  the  center  window 
of  Figure  V-7,  the  side-bands  off  the  main  peak  can  be  seen. 
Due  to  the  power  in  the  side-bands,  limit-cycle  behavior  is 
present  in  the  total  power  P(n),  and  can  be  seen  in  the 
output  power  shown  in  the  lower  right  hand  window.  The 
output  power  oscillates  between  two  power  levels  as  the  FEL 
oscillates  between  two  frequencies.  The  trapped  particle 
instability  does  not  always  cause  limit  cycle  behavior.  In 
Figure  V-8,  the  power  exhibits  chaotic  output. 
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Figure  V-8  Chaotic  Power  Output 

5.    Limit-Cycle  Behavior  in  the  Negative  Tapered  Case 

Figure  V-9  is  the  de synchronism  curve  P(d)  for  the  case 

of  negative  taper,  5=-5.4tc.  For  this  taper  8=-5.4ft,  limit 

cycle  behavior  never  occurred,  or  was  so  small  that  it  was 

not  noticeable . 
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Figure  V-9  P  versus  d  8= -5. 4 

In  Figure  V-10,  at  the  specific  desynchronism  d=0.005 

and  negative  taper  8=-5.4  we  see  there  are  no  trapped 

particles.  Above  a  certain  value  for  d,  d=0.004,  the  trapped 
particle  instabilities  do  not  occur.  This  results  in  smooth 
power  output  with  no  side  bands  around  the  optical 
wavelength  as  seen  in  the  center  window.  The  overall  effect 
is  a  lower  power  output  overall  but  provides  for  cleaner 
output  power.  This  may  be  important  for  some  applications 
where  a  narrow  spectrum  is  important.  For  other  applications 
more  power  may  be  important . 
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****      pel  Pulse  Evolution   **** 
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Figure  V-10  No  Trapped  Particles 

In  Figure  V-ll,  the  effects  of  the  trapped  particle 

instability  can  be  seen  in  the  left  and  center  windows. 
There  is  no  limit  cycle  behavior  as  is  seen  from  the  steady 
power  ?  (n)  window  in  the  bottom  right-hand  window.  In  every 
caper  studied,  there  were  values  of  d  for  which  the  trapped 
particle  instability  occurred  but  no  limit  cycles. 
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****   FEL  Pulse  Evolution  **** 
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Figure  V-ll  Trapped  Particles  with  Stable  Power 

D.     CONCLUSION 

The  data  collected  from  simulations  of  the  Darmstadt 
FEL  describe  how  decreasing  electron  beam  quality  would 
decrease  that  observed  gain.  Also  simulations  showed  that 
there  are  regions  of  operation  where  limit -cycle  behavior 
and  chaos  may  be  expected  to  occur  in  the  experiment.  If  was 
also  observed  that  the  -rapped-particle  instability  may 
occur  without  the  presence  of  limit-cycles. 
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